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ABSTRACT 



Hirota's method is used to construct multi-soliton and plane-wave solutions for affine Toda 
field theories with imaginary coupling. 



1 Introduction 



There has been a revived interest in Toda field theories [1-12] recently, chiefiy because they 
represent a large class of integrable two-dimensional models which include not only conformal 
field theories (CFT's) but also massive deformations away from conformality. Toda field 
theories (TFT's) have a Lagrangian realization (and hence provide one for the CFT's which 
can be written as TFT's) and can be classified according to their underlying Kac-Moody 
algebra g. If (7 is a finite-dimensional, semi-simple Lie algebra, the theory is a CFT. If 
g is an affine Kac-Moody loop algebra, the resulting affine Toda field theory is massive, 
but still integrable. Indeed, these theories provide explicit realizations of the phenomenon 



investigated by Zamolodchikov |]13[ concerning integrable massive deformations of CFT's. 
Finally, it has been observed that if one has the full affine Kac-Moody algebra (i.e., 
with central extension and, in fact, two more fields compared with the usual TFT), then the 
theory is conformal again and is called a conformal affine Toda field theory. The affine TFT 
is then considered a "gauge-fixed" version of the conformal affine theory. 

In this paper we are concerned with finding exact solutions to classical affine TFT's. 



To this end we use the method introduced by Hirota [T^ and applied to TFT's in the 
stimulating work of Hollowood |l[]. If the coupling in the theory is real, then the Toda fields 
can be consistently taken as real fields. The theories are unitary and the S-matrices have been 



solved [|1^]. Since the interaction potential has a unique minimum, there are no solitons. But 
these real affine TFT's are not those related to unitary CFT's. Curiously enough, it turns 
out that to have the relationship to CFT's mentioned above, one must take the coupling 
constant to be purely imaginary [|l|. But then one has the additional feature that the Toda 
potential is periodic, i.e., invariant under translations on the weight lattice of g. Hence, the 
vacua are infinitely degenerate and there is a rich spectrum of topological solitons. But the 
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Lagrangian is complex and the Hamiltonian is non- Hermit ian; so the theory would appear 
to be non-unitary. Yet the solitons have real masses, and it has been suggested that for 
some region of the imaginary coupling, there may exist a truncation which would render the 
theories unitary 0. 

Hirota's method has been successfully applied to the affine Toda theories with purely 
imaginary coupling. Single-soliton solutions were obtained in [|l], |^, ^ and multi-soliton 
solutions for A^^^ theories in In this paper, we will construct more multi-soliton solutions 
to the complex affine Toda theories. 

In §2, we will review briefly the affine Toda theories and Hirota's method. Since the 
solutions for non-simply-laced Toda models can be derived from the simply-laced ones, we 
will mostly work on the simply-laced cases with some comments on non-simply-laced cases. 
In §3, we construct multi-soliton solutions for the Ar loop algebras, in §4, we deal with so(8), 
and in §5, D2r, We summarize our results and conclusions in §6. 

2 Affine Toda Models and Hirota's Method 

The Lagrangian density of affine Toda field theory can be written in the form 

£ = lid.^W^) - ^^E^,(e^"-^ - 1). (2.1) 

The field y^{x,t) is an r-dimensional vector, r being the rank of some finite dimensional 
semi-simple Lie algebra g. The a/s, for j = l,...,n are the simple roots of g; ip is the 
highest root: ip = Z]j=i^i«j; where the n/s are positive integers, and no = 1. ao = —i^ is 
the extended root of g, the affine extension of g. P is the coupling constant, m is the mass 
parameter. Note, in general repeated indices are not summed over in this paper. 
The equations of motion are 



Writing out in component form, = I]5=i(2aj/a|) f-', we have 

^ = -^l, exp(/? K,kv') - exp(/5 E K^kV') , (2.3) 

P V fc=0 fc=0 / 

where i^^-fc = 2a j ■ ak/a\ is the the Cartan matrix of and //s are integers related to n^'s 
by /j = {a'^/ip'^) rij, for j = 0, 1, ■ ■ ■, r. As noted in the introduction, we are interested in 
the case where the couphng constant is purely imaginary. 

We now use Hirota's method to solve the equations of motion. One introduces a nonlinear 
transformation 

= -— E In r^-, or ip^ = -7-^ (in Tj - Ij In Tq) , (2.4) 

where the new variables are the tau-functions for the theories. The equations of motion 
( p.2|) then decouple to a collection of equations in Hirota's form, 

^ ( - Dl) r, . r, = /, (n n""^' - 1) (2-5) 
Here and Dt are Hirota derivatives, defined by 

We expand the r-functions in a formal power series in an arbitrary parameter e as 

^. _ .nAn) _ 1 j_ .^(1) j_ .2^(2) 



X— a;' 



= ^ eV(") = 1 + er]^) + e'r^ + ■■■. (2.6) 

n=0 

Using the bilinearity of Hirota's operator, we can rewrite ( |2.5| ) as 

^ E E ^"(A^ - ^^"^ ■ r^'' = lJf[ - 1) rf- (2.7) 

^''t n a+b=n \k=0 1 
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Expanding the r.h.s. also as power series in e, and equating corresponding coefficients, the 
resulting equation can be written in the form: 

E L,,rt^ + - dl) rf = uf, (2.8) 

k ^ 

where Ljk = IjKjk, and f/j"'* collects the rest of the n^^ order terms which depend on r*^^\ 
• • -, r*-""^-'. The explicit form of Uj"'^ depends on the Cartan matrix K. Now the nonlinear 
equations of motion are simplified to a set of linear partial differential equations, which can 
be solved iteratively. Furthermore, the ffist order equation is homogeneous, Uj^^ = 0. 

For a single soliton traveling at velocity v, tj should depend on x and t only through the 
combination z = x — vt, rj{x,t) = Tj{x — vt) = Tj{z). Here equation ( p.8|) simplifies to a set 
of ordinary differential equations: 

E - '-^^2-t\^) = (2.9) 

with the first order equation given by: 

(l-t;^)-^^ )(z) =m2^L,,r«(z). (2.10) 



Equation (|2.1(J|) can be solved by diagonalizing the matrix Lj^. There are two linearly 



independent solutions for any given eigenvector A (with corresponding eigenvalue A). For 
the null eigenvector Aq = (/q, ■ ■ i ^r), with eigenvalue Aq = 0, the two linearly independent 
solutions are: 

T^'\z) = {CI+C2z)lj, (2.11) 

which correspond to vacuum solutions. For A with non-zero eigenvalue A, the two linearly 
independent solutions are: 

rf^(;2) = A,e^^+^, (2.12) 

where ^ is an integration constant, and 7 satisfies (1 — f ^)7^ = m^A. Since 7 can be either 
positive or negative, ( p.ll| ) and (|2.12|) give all the linearly independent solutions to (|2.10|) . 



General solutions to ( f^.lO| ) are superposition of the solutions given in ( ^.121 ), solution to 
can also be superposed with a solution to ( |2.1(J| ), but they will not have properties of 



a single soliton, and will be dealt with later as multi-soliton solutions. Here we will restrict 
our solutions to the single-soliton form: 

= 5fe^^, (2.13) 

where F = 7(x — vt) + ^, and is independent of x and t. Now ( |2.9|) can be simplified to 
a set of linear algebraic equations: 

E (l,, - n^A5,,) 5f = (2.14) 

k 

where V^''"^ depends only on (5*^^^ ■ ■ ■, 6'^'^'^\ In particular, V^^^ = 0, so 5*^^'' is an eigenvector 
of L with eigenvalue A. S^'^\ S^^\ ■ ■ •, can be solved from (|2.14|) uniquely, except for su{6p) 



and sp{Qp), where additional degeneracy exists. In it is shown that 6j^^ vanishes if n > kIj 
for some integer k, which is less than or equal to the degeneracy of the eigenvalue A, except 
for the special cases of su{6p) and sp{6p) noted above. 

Since the equation ( p.8| ) for rj^-* is linear and homogeneous, one can superpose the single- 
soliton solutions to get general N-soliton solutions, at least up to the first order: 

=ei4!K^ + --- + 6^<5gy-, (2.15) 

where for a = 1, • • ■, N are some arbitrary expansion parameters, S^^j is the first order 
coefficient of the a^^ soliton and is an eigenvector of the L matrix, as discussed above. For 
a general N-soliton solution, we use the following expansion for tj: 

'^j = H ■ ■ ■ e^"^ 5j(mi, ma, • • • , thn) exp (miFi + msFg H h mArFAr) , (2.16) 

{rria} 



where 5^(0, 0, ■ ■ ■ , 0) = 1, and from (|2J[5| ) 



5,(0,---,0,m, = l,0,---,0) = 52. (2.17) 
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From now on, we will let = 1. 
Let us define 



/(mi,m2, ■ ■ ■ ,mAr) = ^ (("^i7i H V mN'^Nf - {mi'^iVi ^ VniN'^NVNf) 



= ^mlXs, +2 Y,'m'amb^JXsa>^st cosh Oab, (2.18) 

where 9ab is the rapidity difference of a^^ and b^^ soliton, and A*^"-* is the eigenvalue for the 
a*^ soliton, we have assumed that all the 7's are positive just for definiteness. Using the the 
following property of Hirota's derivative operators: 

(D^ - Dl) e"l"+^l* . e"2-+/32* = ((/^^ _ _ (^^ _ ^^^^'2Ya,+a^)x+{Pl+P2)t ^ ^2.19) 

we can write: 
1 



2m^ 



{Di - Di) r, ■ r, + /,r/ = (2.20) 
- E((/(^1' ■ ■ ■ , - 2/,)5,(fci, ■ ■ ■ , fc^) + A,(fci, ■ ■ ■ , fc^))e*^iri+-+'=-r-, 

where 

Aj(/ci, • • ■ , /cat) = ^ Q/(A;i - 2mi, • • -, /cat - 2mAr) - /jV 

0^(mi ,---,m]v)-<(fci,---,fc]v) 

(5j(mi, ■ ■ ■ , m]s[)5j{ki — mi, ■ ■ ■ ^k^ — nif^). 

The partial ordering (mi,---,m7v) ^ (fci, ■ ■ ■ , /ca?) is defined as < fca, for all a's, and 
S'^a < ^a- Also we can write: 

n ^k""'" = E (- E ^.A(^i, ■ ■ ■ , A^iv) + B,{k^, ■ ■ ■ , fc^))e'=^r^+-+'=-r-, (2.21) 

where Bj{ki, ■ ■ ■ , k^) is a function of the Sj{mi, ■ ■ ■ , m7v)'s for (mi, ■ ■ ■ , m^r) -< (fci, ■ ■ ■ , kjy). 
The explicit form of Bj{ki, ■ ■ ■ , /ctv) will depend on the Cartan matrix K. Using ( p.2(]| ) and 
( p.21| ), the nonlinear partial differential equation ( |2.5| ) can now be simplified to a set of linear 
algebraic equations: 

E(-^ifc - f{ki, ■ ■ ■ , kN)Sjk)Sk{ki, ■ ■ ■ , /cAf) = Vj{ki, ■ ■ ■ , /cat), (2.22) 



where Vj{ki, ■ ■ ■ , k]\f) = Aj{ki, ■ ■ ■ , + ljBj{ki, ■ ■ ■ , k^)- Since the r.h.s. of ( [^.221 ) only 
depends on the preceding 5j's, as discussed after ( ^.211 ), the above equation can be solved 
iteratively. 

First we decompose the right hand side of (p.22|) into linear combinations of the eigen- 
vectors, 

Vj{ki, ■ ■ ■ , fc^) = ^ Cs{ki, kN)Asj, (2.23) 

s 

where A^, for s = 0, 1, ■ • ■, r, is the complete set of eigenvectors of the matrix L. Since L 
is not necessarily symmetric, the eigenvectors are not necessarily orthogonal to each other 
with respect to the usual inner product. But for simply-laced loop algebras, eigenvectors of 
the matrix L corresponding to distinct eigenvalues are orthogonal to each other with respect 
to the diagonal metric (l//o, l/^i, ■ ' ' ; This observation will make the decomposition 

( p.23| ) easier. Now from equations (|2.22|) and ( p.23|) one can see that 



s n 1 \ sr^ Cs{ki, ■ ■ ■ , k]\f)As j /o o/i\ 



Because of the iterative nature of ( p. 221 ), and from the initial choice ( ^.17| ), one can show 
that 

6j{0, ■ ■ ■ , 0, m^, 0, ■ ■ ■ , 0) = 6a,jima), 

6j{0,- ■ ■ ,0,ma,0, - ■ ■ ,mb, - ■ ■ ,0) = Sab,j{ma,mb), (2.25) 



where the Sab,j{ma,mf,) is a coefficient from the pair scattering solution of the a^^ and b^^ 
soliton, etc. This means that the coefficients of (A^ — l)-soliton solutions can be used in the 
N-soliton solutions, this property makes solving general N-soliton solutions tractable. 

We conclude this section with a number of observations and remarks, the first of which 
concerns integrability. Due to it, the perturbative expansion terminates at finite order. By 
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counting the powers in e, it is shown in |^ that, for a single-sohton solution, the highest 
nonvanishing order in Tj is nlj for some constant integer k. Assuming that the expansion 
( p.l6| ) gives exact solutions, then by repeating the power counting for each e^, one can show 
that: 

6j{mi,- ■ ■ ,mN) = 0, if (mi, ■ ■ ■ , mA^) 7< (k^/j, ■ ■ ■ , k^/j). (2.26) 

For all the double-soliton solutions we obtained, we can explicitly verify that k,'^ = Ka for 
a = 1,2, where Ka is the k value for the a^^ soliton. For a general N-soliton solution it is 
also true that k'^ = Ka for a = 1,2, ■ ■ ■ , N, since otherwise the resulting expression would 
not reduce to N single-soliton solutions in the asymptotic limit. Let Imax = max(/i, ■ ■ ■ , 1^), 
from the above discussion and (|2.22|) we get: 



V,(mi,---,mjv) = 0, if {mi, ■ ■ ■ ,171^) 7^ Imax {i^i, ■ ■ ■ , i^n)- (2.27) 

We consider a solution found if all nonvanishing 6j^s are found. 

If the Dynkin diagram of the loop algebra gi has non-trivial automorphisms, then gi 
can be twisted to form a new loop algebra §2- The equations of motion of the Toda theory 
associated to gi will also have related symmetries, and thus can be reduced to the equations 
of motion of the Toda theory associated to g2 by identifying certain fields. Also, solutions to 
the gi theory may be reduced ("folded") to solutions to the §2 theory [T6|. Furthermore, 



every solution to the §2 theory can be unfolded to a solution to the gi theory, so we only 
need to deal with simply-laced loop algebras. 

Finally, we want to point out how plane-wave solutions come out of the formalism. If we 
replace the single-soliton ansatz (|2.13|) by the plane-wave ansatz 

^^Wgin*^ (2.28) 

where $ = {kx — uot) + ^, and u"^ — k"^ = m?\, most of the discussion on the solitons is still 
valid for the plane waves. Even though we will only write down explicitly the soliton — soliton 



scattering solutions, soliton — plane-wave and wave-wave scattering solutions can be easily 
derived by changing some parameters from real to purely imaginary. 

3 Ar theory 

The Dynkin diagram for A^^^ is shown in Figiyyre 3.1. 

o 

Figure 3.1: Affine Dynkin diagram for A^. 

Since Ij = 1, for all j's, so Lj^ = Kjk, and Hirota's equation ( [^.5D reads 

^(A' - Dl) r, ■ r, = (r,+xr,_i - r/) , for j = 0, 1, 2, ■ • • , r. (3.1) 

Due to the periodicity of the extended Dynkin diagram, it is implied that Tj+^+i = Tj- The 
explicit expression for Bj{ki, ■ ■ ■ , k^) is given below: 

Bj{ki,- ■ ■ ,kN) = X] Sj+i{mi,- ■ ■ ,mN)Sj^i{ki - mi,- ■ ■ ,kN - rriN)- (3.2) 

0^(mi ,---,m]v)-<(fei,---,fc]v) 

The matrix L has eigenvalues and eigenvectors 

K = 4smH,, A,= (l, e^'^% e^''^%---, e^-''^), for s=l, 2, ■ ■ -, r, (3.3) 

where -^s = sn/h, and h = r + 1 is the Coxter number. 

For any given integer 1 < s < r, there is a class of single-soliton solutions (k = 1) 
given in [|l|, with the only non-vanishing coefficients S^^^ = Ag. The topological charge of 
the soliton depends on the value of the parameter ^. Because of degeneracy in eigenvalues, 
\s = Xh-s, there also exist k = 2 solitons 0: 

6j^^ = exp {2msj/h) + exp {—2Txisj/h) , = y^y^ cos^ d^. (3.4) 
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ar-2 Otr-1 



Here and y~ are independent variables, if we take either of and y^ be zero in the above 
equations, we will get back to k = 1 solutions. In the following we will construct multi- 
soliton solutions consisting of solitons with k = 2, since if there are some k = 1 solitons 
involved, we only need to take some of the y's to be zero. 

First we solve for a two-soliton solution. Let A*^^) = Ap, A*^^-* = Ag, p, q are integers with 
1 < P, q < r, and Ki = k,2 = 2, from ( p.l5| ), ( p.25|) and 



0) = yt exp {2nipj/h) + y^ exp {-2nipj/h) , 
^ji^A) = yt exp {2mqj / h) + y2 exp {-2mqj / h) , (3.5) 

^jC^^O) = ytVi cos'^^p, S j {0, 2) = y^^y 2 cos'^'& q. 
By straightforward calculations, we get: 

'^j(l) 1) = ytvt e^^"'' exp (27ri(p + q)j/h) + y^y^ e>« exp (-27ri(p + q)j/h) 

+ yt y2 e'^^' exp (27ri(p - q)j/ h) + y^y^ e^""^ exp (-27ri(p - q)j/ h) , 
5,(2,1) = 5,(2, 0)5,(0, l)exp(7++ + 7pV)' (3-6) 

5,(1,2) = 5,(0, 2)5,(1, 0)exp(7++ + 7+-)' 
5,(2,2) = 5,(2,0)5,(0,2)exp(27++ + 27+-), 

and all higher order terms vanish. Here 

/ ^4- X / X / ^ X cosh ^af, — COS _s, „, 

exp 7++ J = exp 7;,:,J = exp 7+,,, = — rf^ 3.7 

" " " cosh a^ab - cos 

( +- \ I -+ \ ( - \ coshe„fc + cos^9,„+,, 
exp(7.t,.J = exp(7s„,l) = exp(7. .J = ^ , • (3-8) 



cosh^^o?, + cos-i?, 



Constructing an N-soliton solution means finding all nonvanishing 5,- (mi, ■ ■ ■ , mAr). From 
( p.25| ), if any of the vanishes, then the corresponding 5 can be taken from some — 1 



soliton solutions. Since we have assumed k = 2, from (|2.26|) , if any of the nta is greater 



than or equal to 3, then the corresponding 5 shall vanish. So we only need to consider 
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1 < f^a ^ 2. Furthermore from the above two-sohton solutions we see that 5,(2,1) and 
2) are related by renaming the solitons, or shuffling the indices p ^ q. This is also true 

i N-i 

for N-soliton solutions, so we only have to find 5,(2, ■ ■ ■ , 2, 1, ■ ■ ■ , 1), for i = 0, 1, ■ ■ ■, N. 

Let A^"-* = Xs^, with 1 < Sa < r, We list the nontrivial coefflents for N-soliton solutions 
below: 

TV 

5,(T~^) = 5:yr---2/7exp(227r(aiSi + --- + a;vSTv)j-A)exp( J2 ^Zsl)^ (3-9) 

{a} l<b<c<N 
i N-i i N-i i N-i 



5,(2,---,2,l,---,l) = 5,(2,.--,2,0,---,0)5,(0,---,0,l,---,l)nexp (7^,., + 7^".,.,), (3.10) 



b>i 



N 

5,(2~^) = 5S5 ■ ■ ■ 5{,^). n exp (27+,., + 27;,,.,). (3.11) 

a<b 

We will outline our construction and notation in three steps, assuming all the — 1 
soliton solutions are known. 
Step 1: To get 

N 



One can easily see that 5, (1, •••,!) should be N-linear in the ?/'s. Note that ytyi is 
considered to be quadratic while yiy2 is bilinear. So 

N 

5, (T~^) = ^ C(cri, • • • , cr;v)z/r • ■ ■ 1/7, 
W} 

where 0"^ = ±1 and C(cri, ■ ■ ■ , ctat) are independent of x, t and all the y's. ?/°" should be 
understood as y^ or y^ , Summation is over all possible combinations of a's 

If we take y^ = ■ ■ ■ = y]^ = 0, then only one term survives, and all the k's become 1 
instead of 2, so we should get those of HoUowood, [| i.e., 

C(l,---,l)=exp( Yl 7+sJexp(i27r(si + --- + s^)j//i). (3.12) 

l<b<c<si 



^In ^ the N-soliton solutions were written down without explicit proof. 
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Proving the above equation is equivalent to proving the following identity: 
^(/((Ji, 0-2, ■ ■ ■ , ctat) - 4 sin^ (ai^^i H h (Tni^n)) Yl {coshOab - cos {aai!^a - crbA)) = 0, 

{a} a<b 

(3.13) 

where each a takes value 1 or —1, and summation is over all possible combinations of a. 
Since directly proving the above identity is hard, we will construct (|3.12|) without direct 



proof of (|3.13|) , by using ( p.27|) instead. Since we have taken yi = ■ ■ ■ = yjy = 0, all the 



single solitons have k = 1. By ( p. 261) , the only nonvanishing terms in 2, ■ ■ ■ , 2) are 



1, ■ ■ ■ , 1) (5,+i(0, 1, ■ ■ ■ , 1) + 1, ■ ■ ■ , 1) (5,_i(0, 1, ■ ■ ■ , 1) 

-2 cos(27rsi//i) 6j{l, 1, • • ■ , 1) Sj{0, 1, ■ ■ ■ , 1). 



Let 5j{l, 1, ■ ■ ■ , 1) = e'2''^("2+- +«iv)/'^T_,-, from the vanishing condition of V,(l, 2, ■ ■ ■ , 2) we get: 

Tj_i + T,+i - 2 cos{27T si/h)Tj = 0. (3.14) 

The above homogeneous second order difference equation has two linearly independent so- 
lutions: 

From a symmetry argument or the vanishing conditions Vj (2, ■ ■ ■ , 2, 1) = 0, we can see that 
we should chose the "+" sign in the above solution. The constant factor T can be determined 
from the vanishing condition of V^(l, 1, 2, ■ ■ ■ , 2), which gives ( |3.12| ). We could find all the 
C(cri, ■ ■ ■ , o"Ar)'s the same way, thus proving ( p.9| ). 
Step 2: To get ( PTTUQ . 



First let = ■■■ = y^^ = 0, so ni = ■■■ = hii = 2, Kj+i = ■ ■ ■ = = 1. Adding 

i N-i-1 



all the nonvanishing terms in V^ (4, ■ ■ ■ , 4, 2, ■ ■ ■ , 2, 1), from Eq. ( p.27|) , the sum should van- 



ish. So we obtain a homogeneous second order difference equation of 5^(2, ■ ■ ■ , 2, 1, ■ ■ ■ , 1). 
We can write down the vanishing condition of V, (0, ■ ■ ■ , 0, 2, ■ ■ ■ , 2, 1) also, which gives a 
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homogeneous second order difference equation of 5j(0, ■ ■ ■ , 0, 1, ■ ■ ■ , 1). Comparing the two 
equations, and from the fact that they obey the same boundary conditions and have similar 



symmetries, we see that Sj{2, ■ ■ ■ , 2, 1, ■ ■ ■ , 1) should be proportional to 5^(0, ■ ■ ■ , 0, 1, ■ ■ ■ , 1). 
The proportionality constant can be determined from the vanishing condition: 

i N-i-2 

y,(47^,27^,l,l) = 0. 

Second, there should be 2^^* terms, and each can be derived in a similar way. Putting 
them all together, we will get (p.lO|) . 



Step 3: To get ( pT 



From the vanishing condition of V{4, ■ • • , 4, 3) and using factorization of (|3.10| ), we obtain: 



(Ajv - 2)5,(2, ■ ■ ■ , 2) + e±*2-'^/'^5._,(2, ■ ■ ■ , 2) + e^'^"" 5 j+i{2, ■ ■ ■ , 2) = 0. (3.15) 

Adding and subtracting the two equations, we see that 5,(2, ■ ■ ■ , 2) is constant, which con- 
stant can be calculated from the vanishing conditions of y(4, ■ ■ ■ , 4, 3, 3). 

This completes the construction. These k = 2 solitons are not really "single solitons" 
in the strict sense, since they change shape after scattering with other (k = 1) solitons. In 
general a k = 2 soliton has two peaks, and each peak behaves as if it is a free soliton when 
scattering off other solitons. These k = 2 solitons should be interpreted as compounds of 
two K, = 1 solitons interacting in a peculiar way so that there seems to be no interaction 
among the peaks. Indeed, one can obtain such k = 2 soliton solutions by taking the static 



limit of some two-soliton scattering solutions of k = 1 solitons . 

Let us look into the scattering of two k = 1 solitons in a little more detail. Since the 
only effect after scattering is a shift in the center of mass for each soliton, and the shift does 
not depend directly on the topological charge of each soliton but rather on the topological 
classes they belong to, one may be tempted to believe that the interations among the solitons 
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depend only on the classes. Furthermore, since 7^ in ( |3.7| ) is always negative, it appears 
as if the interations are always attractive. However, on closer inspection, one notices that 
the two solitons could either speed up or slow down as they approach each other. If they 
speed up, eventually the two peaks will merge together and then separate again. If, on the 
other hand, they slow down, as they approach to within a certain distance, they begin to 
exchange shape, topological charge and other attributes, and then move apart again as if 
the two solitons have switched position. Because the two solitons performed this exchange 
over some finite distance, even though the interaction is repulsive, the total shift looks like 
an attractive case. In both cases, there is no reflection even if the masses are equal. Details 
of which pair of solitons attract or repel each other still need to be worked out. Since we can 
also have two solitons stand still as if there is no interaction at all, the interaction cannot be 
described by a simple equivalent potential. This may be because the solitons are extended 
objects. 

If r = 2N — 1, for some integer N = 2, 3, ■ ■ ■, and 

Tj = Tr-j, for j = 1, 2, ■ ■ -, N, (3.16) 

then the set tq, ti, ■ ■ ■ , tn will form a tau-function solution set for sp(2N) affine Toda 
models. It is easy to check that if we set = y~, then the n = 2 soliton given in ( |3.4| ) 
will satisfy the condition ( |3.16| ). These k = 2 solitons of sp(2N) do look and behave like 



single solitons: the two overlapping peaks do not drift apart since they are protected by the 
symmetry ( |3.16| ). From the solutions ( |3.9| ) to ( |3.11| ) one can also easily obtain multi-soliton 



solutions for sp(2N) affine TFTs. 

4 D4 

Dynkin diagram for D4 is shown below: 
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o o 

Figure 5.1: Affine Dynkin diagram for D4. 
lo — h — h — — 1) h — 2. 

Boih, • • • , /c^v) = Bi{ki, • • • , /cjv) = -63(^1, • • • , /cjv) = B^ih, • • • , A;^) = 0, 
Biih, • • • , /cjv) = ^ So{mi, mjv)5i(ni, . . . , njv)53(pi, . . . ,pjv)54(gi, • • • , ^jv), 

where the summation is over Pa + Qa + + n,a = ka- The eigenvalues and eigenvectors of 
the matrix Lj^ are: 



Ai 


= 2; 


Ai = 


(1, -1, 0, 1, -1 


A2 


= 2; 


A2 = 


(1, 1, 0, -1, -1 


A3 


= 2; 


A3 = 


(1, -1, 0, -1, 1 


A4 


= 6; 


A4 = 


(1, 1, -4, 1, 1), 



We also define Ag = (0, 0, 1, 0, 0) for later convenience. 

For A = 6, the eigenvalue is not degenerate, there is only one corresponding single soliton 
solution with k — 1: 

5(1) = A4; 5(2) = Ag. 
However, the eigenvalue A = 2 has degeneracy 3, so one can construct a k — 3 soliton: 

S^^^ = yiAi + ?/2A2 + Z/3A3, 

<^^'^ = ^{ymM + ?/2Z/3Ai + ymM) + {yl + yl + yD^e, 

5(3) = P3(A4-12A6), 5(4) = P4A6, (4.17) 
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where P3 = 7/17/2I/3/27, P4 = + ?/^?/| + ?/2l/3/9- O^^^ can obtain three k = 2 or k = 1 
sohtons if one or two of the y's is taken to be zero. 

We will first study the two-soliton scattering solutions. Denote: 

_cosh0-l _2cosh0-l _2cosh0-V3 _2cosh^ + l 

' ~ cosh^ + 1' ' ~ 2cosh^ + l' 2coshe + y3' 2cosh^-l- ^ ^ 

Since it is straight forward to calculate the double-soliton solutions, we will write down the 
nontrivial coefficients without derivation. 
For A(i) = A(2) = 6: 



(5(1,0) = (5(0,1) = A4, 

2(5cosh^-4) 
^^''')= 3(cosh^ + l) + 

5(l,2) = 5(2,l) = -4CiC2A6 



5(2,0) = (5(0,2) = A6, 
(11 - 14 cosh ^) 
' 3(2cosh^ + 1) 



■A 



4, 



5(2,2) = CfC2^A 



6- 



(4.19) 



For A(^) = 2 with k = 3, and A^^^ = 6 with k = 1: 6{m, 0) will be taken as 6^""^ in (|4T7| ), 
5(0,1)=A4, 5(0,2) = Ae, and 



5(1,1) 


= C3(?/iAi+?/2A3 + ?/3A3), 


5(2,1) 


= ^C|(l/il/2A3 + Z/2I/3A1 + Z/1I/3A2) - 4^3(1/? + yl + y3')A6, 


5(2,2) 


= C|(i/2 + y2 + y2)Ae, 


5(3,1) 


= P3C|A4 + 4P3(9C3' + 9C3-1)A6, 


5(3,2) 


= 16P3CIA6, 


5(4, 1) 


= -4P4C32A6, 5(4,2) = P4C|A6, 


5(6,1) 


= -AP^ClAe, 6(6,2) = P^ClA,. 


For A(i) = A(2) 


= 2 with Ki = 3, K2 = 1: Let 5(m, 0) will be taken as 5*-™'' 



(4.20) 



5(0,l)=Ai, 5(0,2) = Ag, and 

5(1,1) = 2/l(CiC4Ao-^|^|^|^A6)+C2(2/2A3 + y3A2), 
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5(2,1 
5(2,2 
5(3,1 

5(3,2 
5(4,1 
5(4,2 
5(5,1 
5(6,2 



/I 2a 4(2cosh0-3) 
y^y^ is^^^^ + 3(2cosh^ + l) -^^ j + 3^^^^^^^^^ + 
{ylClCl + ylCl + i/32C2^)A6, 

-^C'iC2yiy22/3A6, 

16 Ci 



(4.21) 



27 (1 - 2 cosh 0) (1 + 2 cosh ( 

l{chM + yl) + c'2ylyl)Ae, 
± C2yiyiyi ^ 

81 (l + cosh0)(l + 2cosh^) ^' 
^C^ClylylylA,. 



-ymyd,^, 



From the solutions ( |4.19 ) and ( 4.20 ), with a little analysis, we can see that the solitons 
retain their shape after scattering. However in the solution ( [4.2 1| ), the k = 3 soliton Changes 
shape after scattering, as signified by the fact that the constant yi is treated differently from 
1/2 and y^. The /t = 3 soliton can be thought of as a composite of three k = 1 solitons. 
We will substantiate this claim by writing down the three-soliton scattering solution. Let 
5(1,0,0) = yiAi, 5(0,1,0) = 1/2A2, 5(0,0,1) = ysAa. Then 5(2,0,0) = ylK, 5(0,2,0) = 
l/lAe, 5(0, 0, 2) = yfAe, and 



5(1, 1, 0) = yiy2/^i2A3, 5(1, 0, 1) = 1/1I/3D13A2, 5(0, 1, 1) = y2yzD2zAu 

5(1, 1, 1) = i/ii/22/3^i2/^23^i3Ao - yiy2y^ {2D^2D2zDvi + I6D123) Ag, (4.22) 

5(2, 2, 0) = yhlDl^A,, 5(2, 0, 2) = yhlDl^K, 5(0, 2, 2) = ylylDl,K,, 

5(2, 2, 2) = y2^2'l/3^12^23^13A6, 



where D, 



2 cosh 0, 



ab 



ab 



and D 



1 



2cosh^ab + l " (2cosh^i2 + I)(2cosh023 + l)(2cosh0i3 + 1)' 

if all the three-solitons have the same velocity, then 9i2 = O23 = dis = 0, and we can have 

Fi = r2 = Fa = F. Then the above three-soliton scattering solution will have the single 

soliton form with 5*^*"^ = Z]mi+m2+m3=m "^("^i? "^2, "^3), which gives the solution set ( [4.17| ). 



-. Now 
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Because of the non-trivial symmetries of the Dynkin diagram, the solutions can be folded 
to yield solutions to G2 affine TFT by demanding Ti = r-^ = T4. The k = 3 soliton with 
1/1=1/2 = 1/3 and the A = 6 soliton survive the folding procedure. Protected by the discrete 
symmetry, these will be true single solitons of G2 affine TFT. 

5 D^J Theory for r > 5 

The Dynkin diagram for Dl^^ (r > 5) is shown in Figure 4.1. 

o o 
p — o— — o — 

a2 tta ar-3 

a o 

Figure 4.1: Affine Dynkin diagram for DI}\ 

^0 = h = ^r-i = = I5 = 2 for 2 < j < r — 2. The explicit form of Bj{ki, . . . , k]\f) 
is given below: 

Boih, . . . , A;^) = Bi{ki, ...,kN)= Br-i{ki, . . . , /c^) = Briki, . . . , /cat) = 0, 

B2{ki, . . . , k^) = 5o(pi, ■ ■ ■ ,PAr)5i(gi, ■ ■ ■ ,gAr)(53(mi, • ■ ■ ,mAr), (5.1) 

Pa+qa+ma=ka 

Br-2iki, . . . , kj^) = J2 5r{pi,- ■ ■ ,PN)5r-i{qir ■ ■ ,QN)Sr-3{mi,- ■ ■ ,mN), 

Pa + qa + ma = ka 

and the rest of the Bj{ki, . . . , kj^Ys are the same as in (|3^ ). 
The matrix L has eigenvalues and eigenvectors 

A, = 8 sin^ i^s, for 1 < s < r - 2; 

A, = (1. 1, 2^^, 2^^, . . . , 2 '°'^''- (-ly. i-iy). (5.2) 

\r-X = \r = 1\ A._l = (1,-1, 0,0,..., 0), A, = (0,0,..., 0,-1,1), 

where "dg = siT/2{r — 1) = six/h. 



First we establish some useful identities for later calculations of double-soliton solutions: 

r-2 



^cos(2j-l)^9, 

J=2 



P 



r — 3, if p=0 mod h; 

1 + (— cos "i^p, otherwise 



cos(2j - cos(2j - 1)^, 

^ cos dp cos 



-1 



r — 1 
2 cos^ 



if p ± g = mod h; 



l + otherwise 



(5.3) 



(5.4) 



gcos(2^1).,cos(2.-l)..cos(2.-l).. 
i=2 



COSl^p COsdq COS^r. 



-2 + 



r — 1 



4 coster) COsda COS 



ifp±g±m = mod h; 



otherwise 



^-2 (cos(2j - 3)?9p cos(2j + 1)??, + cos(2j + l)dp cos(2j - 3)?9g) cos(2j - 1)?9, 

^ COSl^p COS-i^g COS 

2 / cos3'i9p cos3'(9p\ (r — 1) cos^ 
cost^p cos-i^p ) co^-dp co^-dq 
iS'dp cos3'i9p\ (r — l)cos/ 



COSl^p COS-i^g COSt^m 

(r - 1) cos2^9p+g 

^— ^ — it m ± 



(p-g) 

2'i9 _ 

— |— , if m ± (p + g) = mO' 



S-i^g cos^^m. 



-(l + (-l)P+'?+'") 



cost?™' 



= mod h; 
id h; 



Sing] 



non-degenerate, k = 1, 



cos3'i9p cos3'i9p\ (r — l)cos/ 

cos -i^p COSi?p j COS'&p COSt^g 

/ ,N^4-„4-™\ /cos3'i9r, cos 3791, \ 
, / . / ^ ^ ^ ^ otherwise. 

^ COS-i/p COS-i/p / 

have been given in 0. For 1 < s < r — 2, the eigenvalues are 



;le-soliton solutions 



5^^) = A„ 5^2) = M = (0, 0, 1, . . . , 1, 0, 0). 



(5.5) 



For the eigenvalue A = 2, which is doubly degenerate, we can have a k = 2 soliton: 

5^^) = ?/iA,_i + yaAr, 

5(2) = ( c, C, d2, 4, ■ ■ ■ , rf.-2, (-I)'-'C, (-l)'-^C ) , 

^(4) = = c2(0, 0,1,1,- ■■,1,0,0), 



(5.6) 



where c=iyj + {~ir-^yl)/2h, and d, = (-l)^ ((/i - 2j + 1)^? + (-l)-(2j - l)yl)/h. If we 
let 1/2 = ±(— l)'"/2?/i, then the above k = 2 soliton will be reduced to a k = 1 soliton. 
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We will next write down some double-soliton solutions without derivation. Let A^^^ — Ap, 
and A^^) = A^, for some integers l<p,Q'<r — 2, then 5^(1, 0) = A^j, 5^(0, 1) = A^j, and 



5,(0,2) = 5,(2,0) 



3-I 

rnfi(i')^ 4- '79. A pnsfi9„ — i9,.'l 



COS -dp COS COS -dp cos ■j?^ 

5.(1,1) = 5._i(i,l) = (-i)^+^5o(M), 

cos'j^p cos'j^g ^ ~'~ cos'j^p cos'j?^ ^ ' ^ 2, . . . , r 2, (5.7) 
5,(2, 1) = ^/,- 5,(1, 0)5o(l,l), 
5,(1,2) = «,• 5,(0, l)5o(l,l), 
5,(2,2)=i.,(5o(l,l))'. 

Let A(2) = Ap for some < p < r - 2, A^^) = 2, and 5,(1,0) = Apj, 5,(0,1) = 
A._i ± {-ly/'^Ar, then: 

5(2, 0) = u, 5(0, 2) = (0, 0, 1, -1, 1, -1, • • • , (-l)^-^ 0, o) , 
5(l,l) = (A,_i±(-l)i+^A,)e^ 

5,(1,2) = (-l)^Ap,,e^ for2<j<r-2, (5.8) 
5,(2, 2) = 5,(0, 2) e^^ 

where e^ = (cosh 9 — sin -dp) / (cosh 9 + sin f9p) . 

Let A(i) = A(2) = 2, 5(1,0) = ^i(A._i + (-1)''/%) and 5(0, 1) = z^iK-i - {-ly/^K), 
then: 

5(2, 0)/zl = 5(0, 2)/zl = (0, 0, 1, -1, 1, -1, • • • , (-l)^-\ 0, 0) 

An 1^ 1 + (-!)'• coshg ^ 

^^^'^^^^^^^ r-1 cosh^ + / ° 

-2 2 cos2 (l + (-l)('-+^+i)) cosh ^ - 2 sin^ , 

^^1 cosh^-2sin^^ + l ^ ^^-^^ 

5(2,2)=5o^(l,l) li. 
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Taking the static limit of the above double-sohton solution, and letting yi = zl + z2 and 
y2 = (— — Z2), it will then have single-soliton form with S^"^^ = Y.mi+m2=m ^("^i; ""^2), 
which gives the k = 2 solution of ( ^.61) . 

We are interested in presenting an N-soliton solution. As in section 3, it will be built on 
some (A^ — l)-soliton solutions, and for simplicity, we assume that all solitons involved are 
of the form (5.51). Here we list all the non-trivial coefficients: 



<5o(l, 1, ■ ■ ■ , 1) = 1, ■ ■ ■ , 1) = n Sabfiil, 1), (5.10) 

a<b 

5,(1, 1, ■ ■ ■ , 1) = 5,_i(l, 1, ■ ■ ■ , 1) = (-l)^^+^^+-+^-5o(l, 1, ■ ■ ■ , 1), (5.11) 

1,1,---,1 =2^ ^ -a -a -a ^- [[ e^^b 5.12 

Pi COS^9,iCOS79,2---COS79,^ ^^j^ 

i 

5,(2~^,1,...,1) = 5,(2,..., 2,0,..., 0)5,(0,..., 0,l,...,l)n'5afe,o(l,l), (5.13) 

b>i 

5,(2,...,2) = 5f 5o'(l,...,l). (5.14) 



In Eqs. ( |5.12| ), ( |5.13| ) and ( [5.14[ ), j runs from 2 to r-2. We sketch our construction below. 



One can easily write down the second-order difference equation for 5,(2, ■ ■ ■ , 2) from the 
vanishing conditions of Vj (4, ■ ■ ■ , 4): 

5|(2, ■ ■ ■ , 2) - 5,_i(2, ■ ■ ■ , 2)5,+i(2, ■ • • , 2) = 0, for j = 3, 4, ■ ■ -, r-3. (5.15) 

We need two initial conditions to solve (^.15 ): 



5o(2,---,2) =0, =^ 52(2,---,2) = -\/o(2,---,2) =52(1,---,1) 
1^2(4,- ■■,4) = 0, =^ 53(2,---,2) = 52(2,---,2)= 52(l,...,l). 



From these two conditions and ( |5.15| ) we see that ( |5.14| ) is valid 



One can use the vanishing conditions of Vj (4, ■ ■ ■ , 4, 1, ■ ■ ■ , 1) or \/j (4, ■ ■ ■ , 4, 3, ■ ■ ■ , 3) to 
write down a second-order difference equation for 5, (2, ■ ■ ■ , 2, 1, ■ ■ ■ , 1), for 3 < j < r — 3. 
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The two initial conditions are: 

<5o(2,---,2,l,---,l) =0, =^ 52(2,- ••, 2,1- ■■,1) = -K)(2,---, 2,1,- ••,!), 
F2(4,---,4,l,---,l)=0. 

Factoring out (5j(2, . . . , 2, 0, . . . , 0) na<j< b ^ab,o{^, 1) from the equations and the conditions, 



and comparing them to V^(0, ■ ■ ■ , 0, 1, ■ ■ ■ , 1), we see that (|5.13| ) is vahd. 

One can write down the second-order hnear difference equations for 5^(1, ■ ■ ■ , 1), for 
3 < J < T — 3. They are in the same form as in the A^^^ case. Choosing the y's in ( |3.9| ) 



to be y^" = e and comparing to ( p.l2| ), one sees that ( |5.12| ) satisfies the difference 



equation. However, this does not prove ( |5.12| ) uniquely; one still has to verify that the 
boudary conditions are indeed satisfied. 

6 Concluding Remarks 

Hirota's method has been used to construct exact solutions for affine Toda theories. We 
have endeavored systematically to explore this iterative method and although we do not 
claim completeness, we have presented new, explicit multi-soliton solutions for all the cases 
discussed, including the detailed procedures for how they were obtained. The double-soliton 
solutions given can be proved explicitly, while general N-soliton solutions are hard to verify 
directly. Scattering solutions of Eq, E-j or E^ affine Toda solitons can also be obtained in this 
way, due to the lengthiness of the expressions, they are not included here. In addition to the 
simply-laced loop algebras which we discussed in detail, we have also given some examples 
of multi-soliton solutions of non-simply-laced ones. These can be solved directly using the 
iterative procedures, or indirectly, as we have done, by twisting the simply-laced solutions. 

For simply-laced affine TFTs, the k = 2 or 3 solitons are actually composites of more 
fundamental k = 1 solitons, and can be obtained by taking the zero-velocity limit of some 
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scattering solutions, even though they have single-sohton form. For non-simply-laced affine 
TFTs, the k = 2 or 3 solitons are just as fundamental as k = 1 solitons. 



After this phase of our work was completed, additional recent work by us |]I9| cast 
considerable doubt that many of the soliton solutions actually have real energy density, 
contrary to the hopes expressed by others While not negating the interest of the 

multi-soliton solutions presented here vis-a-vis the knowledge they give us about the detailed 
classical spectrum of the ATFTs per se, nevertheless, these problems about real energy 
densities may cause difficulties in considering the the physical applications of some of the 
theories. 

The single-soliton solutions have a continuous parameter ^, which is related to the shape 
and topological charge of the soliton. For most choices of .^-value, the energy density is a 
complex function of x and t. By demanding the energy density of a single-soliton be real, 
then ^ can only take discrete values. However, the reality is not always preserved in the 



scattering. All this will be explored further in a future publication [|T9 |. 



After our work was completed, we received the paper of Olive et. al [ly], in which 
a different method, based on Vertex operators, for finding soliton solutions for TFT's is 
developed. It will be interesting to compare explicit results found with their method to our 
results. 
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